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Abstract 



> 



Modelling the free fall (and radiative phenomenology) of a massive particle, charged 

or not, in a static and spherically symmetric black hole is a classic, good relativistic dare 

■ that produced a remarkable series of papers, mainly in the seventies of the past century. 

| Some formal topics about the mathematical machinery required to perform the task are 

unfortunately still not very clear; however, with the help of modern computer algebra 

' techniques, some results can at least be tested and corrected. 

O ' 

1 Introduction 

o 
a -1 

When a massive particle falls in a black hole, it induces changes in the spacetime metric that, 
^jfj| although considered as an "ordinary" perturbation, can produce a phenomenology hardly 

comparable to what it would be in a flat background. This fact, even if dealing with ele- 
mentary concepts of General Relativity, found its first serious applications in a celebrated 
paper of 1957 by T. Regge and J. A. Wheeler [Q that somehow founded the theory of the 
stability of the Schwarzschild black hole. However, due to hard complications in the cal- 
culations, this investigation came to a first end, after a number of intermediate stages (in 
particular [21 El lU 03 ) > oni y m 1970-1974 with the work of F. J. Zerilli, who, in collaboration 
with Wheeler and R. Ruffini Ej , attacked the problem of considering the perturbation of 
the e.m. and gravitational fields produced in the falling in greater detail. In jH], in particular, 
a realistic treatment is made, considering the first-order stress-energy tensor contributions 
of the perturbing particle, but the ideas lying at the base of the angular decoupling of the 
equations still remain cryptical in some of their aspects; it is now possible, also by the help of 
quite popular computer algebra means (Mathematica by S. Wolfram and the tensor calculus- 
dedicated package MathTensor by L. Parker and S. Christensen), to account for some hidden 
aspects of the logic path followed by the authors and fix some errors. 
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2 Formulating the problem 



The perturbation analysis is performed by writing the Einstein's system (in geometric units: 
G = c = 1) as 



where h is the perturbation tensor accounting for the presence of a massive particle, eventually 
endowed with an electric charge, that acts directly on the spacetime geometry in addition to 
the Schwarzschild metric tensor: g^ v = g^ u + where corresponds to the line element 



with v(r) = — A(r) = ln(l — 2m/r), m being the mass of the BH. 

Historically speaking, almost all the ideas developed to perform a pure metric perturbation 
analysis (i. e. with T^ u = 0) in the original paper of 1957 were adopted by the successive works, 
in particular the "harmonic syntax" of the angular terms and the distinction between two 
different kinds of perturbations, belonging to different choices among the angular operators 
that generate the multipole expansion itself. 

3 Nature and aims of the angular bases 

To write the equations in a smart basis, accounting for the angular properties of a static, 
spherically symmetric problem (by this fact allowing the separation of the radial parts), a set 
of harmonic "objects" is originally introduced by Regge- Wheeler, by a mechanism of parity 
splitting and repeated derivations of the scalar harmonics; those objects are seen to form a 
useful basis for rank-2 tensors in the Euclidean 3-dim. space and can successively be modified 
to form a basis of the Minkowski spacetime. They are called "tensor multipole" ( TM) basis. 
Later, Frank J. Zerilli modified the original choice in two of the ten elements and noticed that 
the way those tensors are constructed must be consistent with the procedure adopted by Jon 
Mathews j2], consisting in a chain of external products of elements, starting from the spherical 
vector basis, aimed to the forming of a set which, in its spatial Euclidean version, transforms 
under the irreducible representations of SO (3). This set is called "tensor harmonic" (TH) 
basis. 

The procedure to obtain the TM basis resides on the combined action of the time and radial 
projection, Pt, P r , the covariant derivative V and orbital angular momentum L = — ir x V 
operators on the scalar harmonics, that turn out to be linear combinations of the THs. Ex- 
ploring the relationships between bases that originate from external products of unit vectors 
and those derived from operator compositions, one is soon aware of the different charac- 
teristics of those sets concerning parity and with respect to scalar products in £ c (g) £ c and 
in C2{S 2 ) (respectively, the Euclidean complexified tensor product space and the square- 
integrable functions' Hilbert space on the unit 2-sphere). This topic was formerly treated 
in a couple of papers ( |1UI lllj ) that unfortunately contain a certain amount of errors and 
misprints. Some of those results (those, in particular, about the construction of a new TH 
basis of £,? with improved features) are revisited and presented in a correct, implementable 
form in Appendix A. 





ds 2 = -e v{r) dt 2 + e A(r) dr 2 + r 2 {d9 2 + sin 2 (9# 2 ) 
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4 Scalar, vector and tensor harmonics in £ 



To reach a satisfying definition of rank-2 THs, it is natural to hierarchically construct them 
from simpler objects belonging to the same family: scalar and vector harmonics. 
The basis of the Hilbert space of the square-integrable complex functions on the unit sphere 
embedded in £ 3 , whose elements possess the property of being eigenfunctions of L 2 and L z , is 
indicated by the double-indexed function Yjm(9,4>) (the scalar spherical harmonics), and its 
full expression normally adopted (after having made explicit the Legendre polynomials that 
appear in it) is: 



Q\M\ 



d{cos6)\ M \ 



2 J J\ d(cos9) 



°'' (cos 2 9 - 1) J 



(2) 



e iMcj> 



in which the spherical coordinates are, conventionally, 9 6 [0, tt] (the polar angle referred 
to the z-axis) and 4> S [0, 2ir) (the azimuthal angle referred to the x-axis of a rectangular 
cartesian frame). 

In fact, it is of common experience that the practical problem of separating the angular parts 
in a set of equations that must be reduced to a pure radial form can be efficaciously treated 
once a gauge choice is made that transforms those angular parts in algebraic expressions 
depending only on spherical harmonics of same (J, M) and their derivatives (usually up to 
the second order in 9). All combinations of those objects can be rewritten in terms of, at 
most, two different harmonics, in (J, M) and (J — 1, M), making use of well-known recurrent 
relations between Legendre polynomials of different J, leading to the following formulas: 

° J} ' J cot9Y JM - l(J, M) gsc0Y>_ 1iM (3) 



d9 



d 2 Yj 



M - {J 2 cot 2 9 - [J (J + 1) - M 2 } esc 2 9} Y JM + l(J, M) cot 9 esc 9 Fj_i, M (4) 



d9 2 



with l(J, M) = ^ 2J+1 l { /_ 2 - M2) . 

Those two expressions can be joined together, once noticed that = -M 2 Y JM , m a 

second order partial differential equation satisfied by the Yjm- 

d 2 Y JM 2 a 9 2 Yj M +n® Y JM , nv n fK\ 

—Q02- + CSC 9 + COt 9 —90— + J(J+ 1 ) Y JM = . (5) 

This relation will be used later to write in a more compact form the Zerilli's TMs. 

The second step on the way to the THs was made by Blatt and Weisskopf 0, who introduced 
a basis of the space of the complex vector fields on S 2 , which is built as: 

Yj M (9, <t>) = (1, m, l,ra | J, M) Y lm (9, 4>) e„ (6) 
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where ej symbolizes the generic element of the basis of the complexified Euclidean space £% 
composed by the simultaneous eigenvectors of the spin operators S, S z with S = 1: 

ei = — (x + i y)/v2 
e = z 

e^ = (x-iy)/V2 (7) 

and the symbol ( | ) is the bracketed Dirac notation of a Clebsch-Gordan coefficient. 

After this definition, the rank-2 tensor space £| <8> can be provided a proper basis by 

considering orthonormalized external products of the e*; keeping the formalism of [3]: 

l 

t m= (hp, h™>-n\j,m) e m _ M 0e M (8) 

Now, the harmonic tensor spherical basis, spanning the space of the finite-dimensional irre- 
ducible representations of SO '(3), can be defined as: 

3 

y 'jm( q ^)= E (L,M-m,j,m\ J,M) i£ Y IAf _ m (M) (9) 

m=-j 

The orthonormality in C\(S 2 ) of these objects is guaranteed once a scalar product in this 
tensor space is defined as (the overbar meaning complex conjugation): 

/p 2n r 7r 
T:Sdn= / T pa S pa sin 9 dddcf) . (10) 
Jo Jo 

It is worth noting that such a frame is not orthogonal in £% (g) £% endowed with the scalar 
product (T, S) =T : S, for fixed values of J, M. 



5 Tensor multipoles 



It is not difficult to directly write Zerilli's TM covariant basis in Ai 4 once a proper definition 
of L acting on scalar functions is adopted: 



L/ = -irAV/ c = P -irV/;P 
where En V belongs to a rank-2 component of the well-known Levi-Civita's tensor e: 



(11) 
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(12) 



rj pu being the covariant Minkowski metric in spherical coordinates. With this specification, 
the TMs are defined as follows (the symbol "q" means that only the symmetric part of the 
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tensor is considered): 



(13) 



TM\ = 
TM 2 = 
TM 3 = 

TM 4 = 
TM 5 = 

TM 6 = 
TM 7 = 

TM S = 

TM 9 = 
TM W 



(P t oP t )Y JM 

(P r o P r ) Yjm 

V2(iP t o-p r ) Q Yjm 

m(J) (irP t oV) Y JM 
m(J) (r P r o V) Yj M 

m{J) (-iP t o L)q Yjm 
m(J) (P r o L) Y JM 

n(J) [(P r o Dn + r (L o V)q] Yjm tensor magnetic component 

[(L o D() + 3r (P r o V) + r 2 V o V] Y JM 
- ^ [(L o D() - r (P r o V) () - r 2 V o V] Y JM 



scalar components 



vector electric components 



vector magnetic components 



tensor electric components 



where m(J) = yj jj^, n{J) =- y (J „ 1} J(J+1) (J+2) ■ 

The explicit matrix form, with rows and columns labelled as (t, r, 8, 
the help of (JSJ), can be cast as: 
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where " * " denotes a symmetric component of a tensor. 
Under a geometric point of view, the TMs can be classified as: 

• Space longitudinal (i.e. orthogonal to the unit 2-sphere) elements: 
TM 2 , TM 5 , TM 7 (in Zerilli's notation "a", "6", "c"); 

• Space transverse (i.e. tangent to the unit 2-sphere) elements: 
TM 8 ,TM 9 ,TM W (V, "/","<?"); 

• Time addition elements: 

TMi, TM 3 , TM A , TMq ( W/VYVVW')- 

Being the ideas underlying the formation of this set mainly based on the euclidean proper- 
ties of transformation under L and V, orthonormalization problems were expected to arise 
for the components belonging to the — > M.\ extension: in fact, it is straightforward to 
see that the last three cited multipoles (TM 3i ^q) have their norm (defined by (|1U(0 equal to -1. 



5.1 Regge- Wheeler perturbation tensors 

In the TM basis, h m , h e , the well-known Regge- Wheeler perturbation tensors that split 
Q in magnetic and electric parts, totally decoupled from each other, can be written as: 



(14) 



h' 



m(J) r 



[ - h (t,r)TM 6 + ih 1 (t,r)TM 7 



n(J) r- 



h 2 (t,r)TM 8 
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m(J) r 
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n(J) 
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1 



m 2 (J) 
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TM 10 



and, once the celebrated gauge choices (h 2 = for h m and ho = h\ = G = Q for h e ) that 
share the same names are applied, they reduce to: 



(15) 



m(J) r 



[ - h (t, r) TM G + i hi (t, r) TM 7 



h e = e uir) H (t, r) TMi + e x{r) H 2 (t, r) TM 2 + [-% H x (t, r) TM 3 + K(t, r) TM 10 ) 
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making evident their belonging to different classes of transformations of SO(3). 
It is easier to explain their decoupling in those terms than invoking an "opposite parity" 
feature which is hardly recognizable, once we refer to the eigenvalues (if existing) of the 
operator P that acts on a function defined on S 2 as P(f(9, <p)) = f(ir — 9, (ft + 7r): in the case 
of h e , it is manifestly PQa? ) = P(Yjm) = (— 1) J Yjm for any choice of /it, u, while h m has 
no defined parity, since it is immediately seen from © that P(d e Y JM ) = (-1) J+1 Y JM while 
d+YjM = iMYj M so that P^Yjm) = P(Y JM ) = {-l) J Y JM . 



A Harmonics and multipoles in the euclidean case 

In 1976 M. Daumens and P. Minnaert (Universite de Bordeaux- France) produced two papers 
|1U1 II 1 j about the relationship between tensor harmonics and tensor multipoles in both the 
Euclidean and the Minkowski space. Concerning the euclidean 3-dim. space, to stress the 
fact that no conceptual contradictions arise using one basis instead of the other in separating 
the angular dependance of systems like once a unique definition of the quantum index J 
is adopted, they derived an alternative TH basis which had the property of being traceless 
in all its components but one and parity-defined in the traditional sense, contrary to Zerilli's 
TMs, and, contrary to the formerly defined THs, orthogonal also with respect to scalar 
product in (g) Unfortunately, some of the main formulas and relations concerning 
this topic contain errors and/or misprints. The correct form of this harmonic basis for rank-2 
tensors in a Euclidean three-dimensional space is presented here in a fully implementable form, 
illustrating the linear combinations of the THs and TMs previously defined that generate it. 



Called v the element of S 2 whose coordinates are - — ((— 1) L+1 , i, 1) , the basis can be 
defined as: 



J+j 

X J j l M (9,<f>)= r \/3v' e signW (j, \l\, J,-\l\\L, 0) Yj l M (9, <j>) (16) 



L=\J-j\ 



where < j < 2, — j < I < j . 

It is readily seen that those harmonics have the following parity: 

X j j l M {-u) = {sign(/) + [1 - sign 2 (/)] (-1)^} X j j l M (u) , (17) 

u being a generic element of S 2 . Every multipole of order j is found to be a linear combination 
of tensor harmonics of the same order, following this scheme (that corrects ^0] - Table 1): 



(18) 



yOO _ V J 
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\^JM ^JM A JM) — \ Y JM 1 JM 1 JM 11 

(v22 y 21 v 20 y-2-1 v 2 -2\ _ ( v 2 J+2 v 2 J+1 v 2 J v 2 J-l v 2 J-2\ R 
l A JM A JM ^JM A JM J M } I JM Y JM Y JM r JM X J M ^ 
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where A and B are two orthogonal matrices, confirming that the ()16|) structure of a linear 
transformation with Clebsch-Gordan coefficients effectively corresponds to a rotation: 
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Finally, the linear combinations of those THs that generate Zerilli's TMs can be easily found 
by a simple matricial calculus procedure: 

a) called J the jacobian matrix of the transformation from the cartesian rectangular co- 
ordinates (x, y, z) to the polar spherical (9, 4>) on S 2 : 



J 



/ sin 9 cos i 
sin 9 sin< 
\ cos 9 



COS V COS i 

cos 9 sine 
— sin# 



sin 
sin 

cos 4> 
sin 8 



± \ 



(21) 



and A4 the Euclidean 3-dim. covariant metric tensor 

/ 1 \ 



M 



1 







(22) 



\ sin 2 9 J 



the operator linking the spherical representation T s of a tensor to its cartesian one, T c , 
can be expressed, in the matrix algebra terms, as: 



T s = O cs (T c )=J- L T c JM 



(23) 



(3) 

b) called TM- the Euclidean 3-dim. tensors introduced in t 4, counterparts of Zerilli's 
multipoles of the previous section (as matrices, they are obtained by deleting the first 
row and column from the (2, 5, 7, 8, 9, 10)-labelled ones); 
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it is readily found that: 
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